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ABSTRACT 

Intermediate polars (IPs) are cataclysmic variables which contain magnetic white 
dwarfs with a rotational period shorter than the binary orbital period. Evolution- 
ary theory predicts that IPs with long orbital periods evolve through the 2-3 hr pe- 
riod gap, but it is very uncertain what the properties of the resulting objects are. 
Whilst a relatively large number of long-period IPs are known, very few of these 
have short orbital periods. We present phase-resolved spectroscopy and photometry 
of SDSS J233325.92-)-152222.1 and classify it as the IP with the shortest known or- 
bital period (83.12 ± 0.09 min), which contains a white dwarf with a relatively long 
spin period (41.66 ± 0.13 min). We estimate the white dwarf's magnetic moment to be 
/iWD « 2 X 10'^^ Gcm'^, which is not only similar to three of the other four confirmed 
short-period IPs but also to those of many of the long-period IPs. We suggest that 
long-period IPs conserve their magnetic moment as they evolve towards shorter or- 
bital periods. Therefore the dominant population of long-period IPs, which have white 
dwarf spin periods roughly ten times shorter than their orbital periods, will likely end 
up as short-period IPs like SDSSJ2333, with spin periods a large fraction of their 
orbital periods. 

Key words: stars: novae, cataclysmic variables - stars: binaries: close - stars: white 
dwarfs - stars: magnetic fields - stars: individual: SDSS J233325. 92-1-152222.1 



1 INTRODUCTION 

Cataclysmic variables (CVs) are interacting binary stars 
containing a white dwarf primary star and a low-mass sec- 
ondary star in a tight orbit () Warner,, 1995^ In most of these 
systems the secondary star is unevolved and fills its Roche 
lobe, losing material to the white dwarf primary via an 
accretion disc. About a quarter of CVs contain magnetic 
white dwarfs, and these systems are split into two cate- 
gories: the AM Her stars (polars) and the intermediate po- 
lars (IPs). In polars, the strong magnetic field of the white 
dwarf causes its spin period, Pspin, to synchronise to the or- 
bital period, Porb- The evolution of these objects is expected 
to be strongly afi'ected by interaction between the magnetic 
fields of the white dwarf and the l ate-type secondary star 
jWebbink fc Wickramasinghd |2002| ). In IPs it is believed 
that the lower magnetic field of the white dwarf means that 
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its rotation has not become synchronised to the orbital mo- 
tion. The evolution of IPs and in particular their relation 
to polars is still a debated issue l|HeUieij|200TI : [Norton et~all 



|2004 hereafter NWS04). 

While the total number of confirmed IPs has been 
boosted substant ia lly o ve r the past few yea rs (e.g . 



Woudt fc Warneil [20o3: lAraui o-Betancor et al 



Rodriguez-Gil et al.l 



2004 



2005; Rodn'EUOz-Gil ct al.' 
iBonnct- Bidaud ct al. 200i 
Ide Martino et al.l l2006i ) 



Woudt et al 
|2005|^ 



2004 



20031: 



Schlegell 



Gansicke et al.' '200J; 
Norton fc Tanner 20(3; 
de Martino et al.l 120061 ). only a handful of these have 
short orbital periods (i.e., below the 2-3 hr gap which 
is prevalent in the observed population of CVs). Con- 
versely, the known population of polars mos t ly lie s below 
Porb = 4hr (jWeiabink fc Wickramasinghd |2002| ). This 
has led to suggestions that long-period IPs may evolve 
into short-period polars (see iPattersonI [l994h . A detailed 
theoretical study by NWS04 found that in general this 
should be the case. However, those with a white dwarf 
magnetic moment of fiwD ^ 5 x 10^^ G cm^ and secondary 
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stars with weak magnetic fields will remain IPs, because the 
magnetic interaction between the stars is not strong enough 
to synchronise the the rotation of the white dwarf to the 
orbital motion. This picture of the evolution of magnetic 
CVs has not yet been confirmed observationally because 

very few short-period IPs are known. 

We are conducting a research program (Gansickc"2005|) 
to measure the orbital periods of objects spectr o scopi- 
cally identi fied as CVs (|Szkodv et all bool |2003| . |2004 
I2OO5I . I2OO6I ) by the Sloan Digital Sky Survey (SDSS). 
The main motivation of this work lies in the characterisa- 
tion of an homogeneously identified sample of CVs which 
can be used to investigate the properties of the intrinsic 
CV population. The relatively faint limiting magnitude of 
the SDSS spectroscopic observations means that this sur- 
vey is much less biased towards intrinsically brighter ob- 
jects than previous work. A very high proportion of the 
SDSS CVs which we have studied so f ar are faint short- 
perio d systems (ICansicke et al.l l2006l: ISouthworth et all 
|2006| . hereafter Paper I: iLittlefair et all I2OO6I '). In Paper I 
we found that SDSS J023322.61-h005059.5 was a probable 
IP with Porb = 96.08 ± 0.09 min and Pspin ~ 60min. 
Here we present a detailed analysis of another system, 
SDSS J233325. 92-^152222.1 (hereafter SDSSJ2333), which 
we find to be an IP with Porb = 83.12 min and Pspin = 
41.6 6 mii{B. 

ISzkodv et af] (|2005l ) presented the identification spec- 
trum of SDSS J2333 alongside time-resolved photometry and 
spectroscopy. They found modulation at 21 min in the light 
curve and at 82 min in the Ha and H/3 emission line ra- 
dial velocities, and suggested that this object could be an 
IP with an or bital period c lose to the 80 min period miin- 
mum for CVs (i Kni ggd lioO^ . The SDSS spectrum shows the 
Balmer and He I emission lines characteristic of CVs as well 
emission at Hell 4686 A which is a hallmark of magnetic 
systems. Whilst the Balmer and He I lines have the double- 
peaked profile indicative of an accretion disc, the Hell line 
is single-peaked. However, it is difficult to conclude anything 
from this because the line may be contaminated by emission 
from the CIII and NIII Bowen blend. X-ray emission from 
SDSS J2333 was not detected by ROSAT. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 WHT spectroscopy 

Spectroscopic observations were obtained in 2006 September 
using the ISIS double-beam spectrograph on the William 
Herschel Telescope (WHT) at La Palma (Table[TJ. For the 
red arm we used the R316R grating and Marconi CCD 
binned by factors of 2 (spectral) and 4 (spatial), giving a 
wavelength range of 6300-9200 A and a reciprocal disper- 
sion of 1.66 A per binned pixel. For the blue arm we used 
the R600B grating and EEV12 CCD with the same bin- 
ning factors, giving a wavelength coverage of 3640-5270 A 
at 0.88 A per binned pixel. From measurements of the full 
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Figure 2. Greyscale plots of the continuum- normalised and 
phase-binned trailed spectra around the H« and H/3 emission 
lines. Other Balmer lines give similar result, but trails for the 
helium emission lines are too noisy to be useful. 



widths at half maximum (FWHMs) of arc-lamp and night- 
sky spectral emission lines, we find that this gave resolutions 
of 3.5 A (red arm) and 1.8 A (blue arm). 

Data reductio n was undertaken using optimal extrac- 
tion (lHornelll986l ) as implemented in the pamel^EI code 
(|Marshlll98'9l '). The wavelength calibration was interpolated 
from copper-neon and copper-argon arc lamp exposures 
taken every hour. We removed the telluric lines and fiux- 
calibrated the target spectra using observations of G191- 
B2B, treating each night separately. 

A total of 38 spectra were obtained over two nights 
and the average spectrum is shown in Fig.[T] They show 
strong double-peaked Balmer emission and double-peaked 
emission from the other usual suspects. He I, Fell and 
CallK. There is no hint of the Balmer lines of the under- 
lying white dwarf or of Call emission. However, there is 
significant sing/e-peaked emission at Hell 4686 A, which is 
characteristic of systems containing mag netic white dwarfs 
(e.g., DWCnc. lRodn'guez-Gil et al.ll2004h . The trailed spec- 
tra (Fig.[2| show that the emission lines have double-peaked 
emission with a strong S-wave component. The spectrum of 
SDSS J23 33 is reminiscent of those of the short-period IPs 
HTC am (|Kemp et al.ll2002l), DWCnc (iRodn'guez-Gil et all 
1200411 and V1025 Cen (|Bucklev et allll998l ). The equivalent 
widths of Hq, H/3 and Hell 4686 A are 78 ± 4, 56 ± 3 and 
8±1 A, respectively. The ratio ^'^'^Hffl^ = 0.14±0.02 is typ- 
ical for short-period IPs (e.g., DW C ue) and not unusual for 
short -period non-magnetic CVs (e.g. iThorstensen fc FentonI 
I2OO3II . 



2.2 Calar Alto photometry 

Unfiltered differential photometry of SDSS J2333 was ob- 
tained on the night of 2006 September 13 using the 
2.2m telescope at Calar Alto Observatory equipped with 
the CAFOS focal reducer (Table[T|. We used the blue- 
optimised SITe CCD without binning, giving a plate scale 



^ The reduced spectra and photometry presented in this work will 
be made available at the CDS (http : //cdsweb.u-strasbg. f r/| 
and at ,http: //www. astro .keele . ac.uk/~jkt/] 
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Figure 1. Flux-calibrated averaged blue-arm (loft) and red-arm (right) WHT spectra of SDSS J2333. 



Table 1. Log of the observations of SDSS J2333. 



Start date 


Start time 


End time 


Telescope and 


Optical 


Number of 


Exposure 


(UT) 


(UT) 


(UT) 


instrument 


element 


observations 


time (s) 


2006 09 13 


19 57 


23 58 


CA 2.2 CAFOS 


(unfiltered) 


293 


40 


2006 09 25 


00 46 


04 46 


WHT ISIS 


R600B R316R 


22 


600 


2006 09 25 


23 07 


01 39 


WHT ISIS 


R600B R316R 


15 


600 



of 0.53" px~^. The images were reduced and aperture pho- 
tometry was perform ed using the pipeline described by 
iGansicke et alj (|2004h . which apphes bias and fiat-field cor- 
rections within MiD Asfl and uses the SEXTRACTOR package 
11 Bert in fc Arnouts|[T99g ') to obtain aperture photometry for 
all objects in the field of view. The differential magnitudes 
were converted into apparent magnitudes using the SDSS g 
and r photometry of the main comparison star. 



3 DATA ANALYSIS 

3.1 Radial velocity analysis 

We measured radial velocities from the Balmer emission 
lines by c ross-correlation with sing le and double Gaussian 
functions (jSchneider fc Yound[l980l ). For each line we tried 
a range of different widths and separations for the Gaussians 
in order to verify the consistency of our results (see Paper I 
for further details). We found that the best results for period 
determination were obtained using the Ha line and a sin- 
gle Gaussian of FWHM 1000 km s~^, as expected given the 
prominence of the S-wave emission in the spectra (Fig. [5]). 
For a double Gaussian the best parameters were FWHM 
300kms~^ and separation 2000 kms~^. The choice of emis- 
sion line, measurement technique and Gaussian parameters 
does not have a significant effect on the derived period. Ra- 
dial velocities derived for the He I 4471 A line gave similar 
results but with much larger noise as the line is weaker. 

Periodograms were c alculated from the measured Ha 
radial velocities using the IScargld (|l982l ) method, analysis 



|http : //www. eso . org/projects/esomldas/| 



of variance (AoV; ISchwarzenb erg-Cze rnvlll989ll and ort hog- 
onal polynomials (ORT; [sd iwarzcnb erg-Czernvl [19961 ). as 
implemented within the TS^j context in MIDAS. All three 
types of periodogram contain a strong signal close to 83 min, 
accompanied by one-day aliases, and a few small peaks at 
much shorter periods (Fig.|31). As the three periodogram 
techniques agree well we have restricted further analysis to 
the Scargle algorithm, which is the most appropriate for a 
simple sinusoidal variation. We obtained the orbital period 
by fitting a circular spectroscopic orbit (sine wave) to the 
radial velocities (Fig.|3} using the SBOflfl program, which we 
have previously found to give reliable error estim ates for 
the optimised parameters jSouthworth et al]|2005l ). Whilst 
the radial velocities measured using a single Gaussian have 
a larger scatter than those from the double- Gaussian tech- 
nique, the velocity variation in the core of the emission lines 
is much larger. These radial velocities therefore give the best 
results, and we find a period of Porb ~ 83.12 ± 0.09 min. 

To investigate the probability that our chosen alias is 
the correct one we have both fitted sine curves for each of 
the possible alias periods and conducted bootstrapping sim- 
ulations (see Paper I). For the sine curves we would expect 
the lowest errors and residuals, and the highest amplitude, 
for the correct period. The bootstrapping gives results which 
can be directly interpreted as a probability distribution, but 
there will be some bias away from the correct period because 
the process of resampling the observations leads to a loss of 
temporal resolution, degrading the resulting periodograms. 
Thus bootstrapping will underestimate the probability that 



^ ^http: //www.eso . org/projects/esomidas/doc/user/ 98MQV/volb/node220 .html | 

^ Spectroscopic Binary Orbit Program, written by P. B. Etzel, 



jhttp : //mintaka. sdsu. edu/f aculty/etzel/] 
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Table 2. Results of the bootstrapping analysis and the best-fitting circular spectroscopic orbits for the Ho line, both radial velocity 
measurement techniques and for different possible periods. The orbit adopted as a final result is indicated in bold. The results of the 
bootstrapping analysis are given as the percentage of simulations returning a highest peak close to a particular period. 



Radial velocity 


Orbital period 


Velocity amplitude 


Systemic velocity 




Scargle bootstrap 


measurement 


(day) 


(kms-l) 


(kms-i) 


( km s"-') 


probability (%) 


Double Gaussian 


0.051038 ±0.00011 


51.6 ±7.8 


-2.7±5.6 


34 




Double Gaussian 


0.054111 ±0.00012 


56.3 ±6.4 


-0.8 ±4.6 


28 


11.5 


Double Gaussian 


0.057622 ±0.00013 


58.4 ±6.0 


2.6 ±4.3 


26 


85.6 


Double Gaussian 


0.061608 ±0.00016 


57.1 ±7.0 


6.6 ±4.9 


30 


2.9 


Single Gaussian 


0.051149 ±0.00012 


151.8 ± 17.6 


21.3 ± 12.9 


78 




Single Gaussian 


0.054222 ± 0.00008 


168.2 ±11.6 


27.5 ±8.5 


52 


17.4 


Single Gaussian 


0.057724 ±0.00006 


176.3 ±8.8 


38.9±6.4 


39 


82.6 


Single Gaussian 


0.061713 ±0.00009 


175.5 ±12.9 


52.5 ±9.1 


55 






Figure 3. Scargle periodogram of the Ha radial velocities mea- 
sured using a single Gaussian of width 1000 km s~^ (upper panel) 
and the phased velocities with the best-fitting circular spectro- 
scopic orbit (lower panel). On the periodogram the chosen peak 
is indicated with an arrow. 



Figure 4. Calar Alto 2.2m photometry of SDSS J2333. The top 
panel shows the observed light curve, which displays a clear sinu- 
soidal variation. The middle panel contains an ORT periodogram 
and the lower panel depicts the phased light curve with a double- 
sine fit. 



the best period is in fact the correct one. Both the sine curve 
fitting and the bootstrapping simulations clearly favour the 
83.12 min period. 



3.2 Light curve analysis 

The light curve of SDSS J2333 shows clear sinusoidal vari- 
ation with a large amplitude of about 0.2 mag (Fig.|4)|. An 
ORT periodogram of these data has strong power at one half 
and one quarter of the orbital period (Fig.jJ)) but not at the 
orbital period itself. A visual inspection of the light curve 
suggested that it exhibits a double-humped variation so we 
fitted the function 



■ ai sm 



2n{x + (t>i] 



+ a2 sin 



27r(2a; + 02) 



where the optimised parameters are ai and 02 (the ampli- 
tudes of variation), P (the period) and 4>i and 4>2 (the phase 
zeropoints). This process confirmed the double- humped na- 
ture of the variation (Fig.|3} and yielded a period of 41.66 ± 
0.13 min. 

ISzkodv et al.l |200^ found a periodicity of 21 min in 
their light curve of this system, which was obtained with a 
much smaller telescope (1.0 m) than our data. Aside from 
the fact that the double-peaked nature of the variation was 
not apparent in their data, this periodicity is consistent with 
our own measurement and demonstrates that it is stable. 

We have found that SDSSJ2333 has Porb = 83.12 ± 
0.09 min, close to the minimum orbital period for CVs, and 
a light curve with a strong double-humped modulation of 
amplitude 0.2 mag and period 41.66±0.13 min. We interpret 
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the light variation to be due to accretion heated spots on the 
white dwarf primary, which has a largely bipolar magnetic 
field. Therefore we classify SDSS J2333 as an IP with an 
orbital period below the period gap and white dwarf which 
rotates with twice the frequency of the orbital motion. 



4 DISCUSSION 

In TableO we have collected the observed orbital and 
spin periods for the five confirmed IPs with Porb < 2 hr. 
These are plotted in Fig.O Classification as an IP re- 
quires the presence of coherent variation at the whi te dwarf 
spin period over a significant span of time (e.g. iBucklevI 
I2OOOI ). As the photometric observat ions for a n additional 
two objects, RXJ1039. 7-0507 ( Woudt fc Warner 2003 1 and 
SDSS J023322. 61+005059.5 (Paper I), have only a short 
baseline their IP nature needs additional confirmation. The 
white dwarf magnetic moments, fiwD, of the objects in Ta- 
ble|3]have been obtained from NWS04 under the assumption 
that the systems are in rotational equilibrium, either taken 
from their Table 1 or estimated from their Fig. 2. 

SDSS J2333 is joined by the objects DW Cancri and 



V1025 Centauri in having Psj 



0.5Porb and fiwv 



2 X 10 Gem . Assuming a canonical white dwarf radius 
of 10^ m, this gives a field strength of 2 MG for the white 
dwarf. The properties of EXHya are also similar, whereas 
those of the fifth system, HT Cam, are quite different. 

The majority of IPs wit h Porb > 3 hr have Ps pin ~ 
O.lPorb (|Barrett et al.l 1 19881 : ICansicke et~all [iooi) and 
^iWD ~ 2 X lO^^Gcm^ (NWS04). This magnetic moment 
is strikingly similar to those of SDSSJ2333, DWCnc and 
V1025 Cen. This is an empirical indication that long-period 
IPs with Pspin ~ O.lPorb couservc /iwD as they evolve and 
become short-period IPs with Pspin ~ 0.5Porb. This is consis- 
tent with the evo lu tionary picture of magnetic CVs draw n by 
iPattersonI l|l994 ): IWebbink fc Wickramasinghd (|2002l ') and 
NWS04. 

HT Cam is a short-period IP with a lower /iwD of 
2.7 X 10'^^ Gcm"^. In the scenario we have outlined, an ob- 
ject such as this will result from the evolution of systems 
w hich currently resem ble SDSS J223843.84-H010820.7 (Aqr 1 
in I Woudt et~aLll2004l ) for which fiwD = 2.6 x 10^^ Gcm^ has 
been estimated (NWS04). A similar argument can be ad- 
vanced for RXJ1039. 7-0507, which has ^iwD ~ 10^^ which 
is very similar to that of several known IPs (NWS04 ), for 
exai nple IRXS J062518.2+ 733433 (|Staude et al.ll2003D . 

iKing fc Wvmil have found that there is a large 

continuum of equilibrium spin levels for short-period IPs, 
and that an IP with a specific mass ratio, mass transfer rate 
and orbital period can have any of a wide range of values of 
Pspin/Porb, depending on the value of /iwD 

The theoretical calculations of NWS04 suggest that 
long-period IPs will evolve into polars unless fiwD ^ 5 x 
10^^ G cm'^ and the magnetic field of the secondary star is 
weak. In this case they may become EX Hya-like systems as 
the interaction between the magnetic fields is too weak to 
synchronise the rotation of the white dwarf to the orbital 
motion. 

Whilst there are twelve long-period IPs listed in NWS04 
with /iwD similar to that of SDSS J2333, there are only four 
short-period systems with this property. As short-period 
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Figure 5. Comparison of the orbital and spin periods of IPs 
with short orbital periods. Dotted lines indicate where the spin 
period is 0.1, 0.25 and 0.5 times the orbital period. Those sys- 
tems which are not confirmed IPs are plotted with open circles. 
The error bars are smaller than the points for all systems except 
SDSS J023322.61+005059.5. 



CVs are predicte d to be intrinsically far more common tha n 
long-period ones l|de Koollll99'3 : lKolb|[l993l : |Politano|[l996l ). 
this means either that the vast majority of the long-period 
IPs become polars or that the known population of short- 
period IPs is much less complete than for those with longer 
periods. The latter possibility could easily arise as short- 
period IPs are in gener al intrinsically fainter than long- 
period IPs (|Warneij|l99"5l ). 



5 CONCLUSIONS 

SDSS J2333 was identified as a cataclysmic variable from a 
spectrum taken by the SDSS, which shows strong double- 
peaked Balmer emission, double-peaked He I emission and 
the single-peaked Hell emission which is often found in CVs 
containing a magnetic white dwarf. We have measured its 
orbital period to be 83.12 ± 0.09 min from a radial veloc- 
ity analysis of the Balmer emission lines. It is therefore an- 
other SDSS-identified CV with an orbital period close to the 
minimum period for CVs containing unevolved mass donor 
stars. Its light curve shows a strong variability with a double- 
humped nature and a period of 41.66 ± 0.13 min, which is 
precisely half that of its orbital period. We interpret this as 
arising from hot spots on the surface of the white dwarf pri- 
mary component caused by accretion of matter controlled by 
a mostly dipolar magnetic field. SDSS J2333 is therefore the 
shortest-period example of a relatively rare class of short- 
period intermediate polars. An X-ray detection of its spin 
period is highly desirable to further investigate its IP nature. 
The fact that its spin period is precisely half of its orbital 
period suggests the presence of a physical mechanism which 
is maintaining this as an equilibrium state. 

Four out of a total of five confirmed short-period IPs, 
including SDSS J2333, have a spin period of approximately 
half their orbital period. These systems contain a white 
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Table 3. Orbital and spin periods of intermediate polar CVs with orbital periods shorter than the 2-3 hr period gap. The estimated 
magnetic moments of the systems are from NWS04, either taken from their Table 1 or traced from their Fig. 2. 
* The IP nature of RXJ1039-0507 and SP SS J 023322.61+005059.5 has not been confirme d 

References: (1) Flbvmassian et al.l lll998l); (2) iKernp et al.l Jiooi): (3) Ide Martino et al.l ||2005| '): f4 ) iRodn'euez-Gil et al] ||2004|): fS) 
[Patterson et al.l ||2004|'): (61 iBucklev et al.l lll99^: (7) iHellier et al.l [l99a )- fS't lHellier et al.l l|2002h : (9) iMumfordI (|l967l '): (W) IVogt et aU 
II1980I '): (ll) lHellier fc Sproatj 1 I1992I) : (12') IWoudt fc Warned tosl ). 



System Orbital period (min) Spin period fspin Estimated magnetic References 

(min) (min) Poih moment (Gcm^) 



HT Cam 


85.9853 ±0.0014 


8.58430 ±0.00000 


0.100 


0.3x1033 


1, 2, 3 


DW One 


86.1015 ±0.0003 


38.58377 ±0.00006 


0.448 


1.5xl033 


4, 5 


V1025 Cen 


84.62 


35.73 ±0.05 


0.422 


1.8x1033 


6, 7, 8 


EX Hya 


98.256738 ±0.000006 


67.02688 ±0.00001 


0.682 


~5xl033 


9, 10, 11 


SDSS J2333 


83.12 ±0.09 


41.66 ±0.13 


0.501 


~2xl033 


This work 


RXJ1039.7-0507 * 


94.4597 ±0.0001 


24.062 * 


0.255 


0.9x1033 


12 


SDSS J023322.61+005059.5 * 


96.08 ±0.09 


60±5 * 


0.625 




Paper I 



dwarf with a magnetic moment of /iwd ~ 2 x lO"^"^ G cm3 
(corresponding to a field strength of about 2MG). A rela- 
tively large number of long-period (Porb 3hr) IPs have 
spin periods close to a tenth of their orbital periods and 
MWD ~ 2 X 10^^ G cm3 . From this we suggest that the popu- 
lation of long-period IPs with Papin ~ O.lPorb will conserve 
Hwu during their later evolution and become short-period 

IPs with Pspin ~ 0.5Porb- 
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